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Abstract 



We study the problem of minimal resistance for a body moving with 
constant velocity in a rarefied medium of chaotically moving point parti- 
cles, in Euclidean space R d . The particles distribution over velocities is 
radially symmetric. Under some additional assumptions on the distribu- 
tion function, the complete classification of bodies of least resistance is 
made. In the case of three and more dimensions there are two kinds of 
solutions: a body similar to the solution of classical Newton's problem 
and a union of two such bodies "glued together" by rear parts of their 
surfaces. In the two-dimensional case there are solutions of five different 
types: (a) a trapezium; (b) an isosceles triangle; (c) the union of a tri- 
angle and a trapezium with common base; (d) the union of two isosceles 
triangles with common base; (e) the union of two triangles and a trapez- 
ium. The cases (a)-(d) are realized for any distribution of particles over 
velocities, and the case (e) is only realized for some distributions. Two 
limit cases are considered, where the average velocity of particles is big 
and where it is small as compared to the body's velocity. Finally, using 
the obtained analytical results, we study numerically a particular case: 
the problem of body's motion in a rarefied homogeneous monatomic ideal 
gas of positive temperature in R 2 and in M 3 . 

1 Introduction 

In 1686, in his Principia [9], I. Newton considered the problem of body's motion 
in a homogeneous medium of point particles. He assumed that collisions of 
the particles with the body are absolutely elastic, the medium is very rare, so 
that the particles do not mutually interact, and that initially the particles are 
immovable, i.e., thermal motion of particles is not taken into account. These 
assumptions are not satisfied in the ordinary conditions "on earth" , but can be 
approximately valid when considering motion of high-speed and high-altitude 
flying vehicles such as missiles and artificial satellites. 



Newton considered the problem of finding the shape of body minimizing 
resistance of the medium to the body's motion. He solved this problem in the 
class of convex axially symmetric bodies with the axis parallel to the body's 
velocity, of fixed length along this axis and with fixed projection on a plane 
orthogonal to the axis. Due to convexity of the body, each particle hits the 
body at most once, and this fact allows one to write down an explicit analytical 
formula for resistance. The body of least resistance found by Newton can be 
described as follows: the rear part of its surface is a flat disk, which is at the 
same time the maximal cross section of the body by a plane orthogonal to the 
symmetry axis. The front part of the surface is composed of a smaller disk in 
the middle and of a strictly convex lateral surface. 

Let us also mention the two-dimensional analogue of Newton's problem. 
Consider a class of convex figures in R 2 that are symmetric with respect to 
some straight line and have fixed length along this line and fixed width; it is 
required to find the figure from this class such that resistance to the motion 
of the figure along this line is minimal. If the length does not exceed the half- 
width, the solution is a trapezium with the angle 45° at the base; elsewhere, the 
solution is an isosceles triangle. 

Since the early 1990th the interest to Newton's problem revived. In partic- 
ular, there were obtained interesting results related to minimization problems 
in wider classes of bodies obtained by withdrawing or relaxing the conditions 
initially imposed by Newton: axial symmetry [7], |S] and convexity [3], 

9, ma, nu 

On the other hand, the assumptions of absolutely elastic collisions and of 
absence of thermal motion in the medium are, at the best, true only approxi- 
mately. (Note that "absence of thermal motion" means that the mean velocity 
of thermal motion of particles is negligible as compared to the body's velocity.) 
In [5], the problem was studied under the more realistic hypothesis of presence 
of friction at the moment of collision (so that collisions are not absolutely elas- 
tic). In the present paper, we address the minimization problem in a medium 
with thermal noise of particles. 

A convex and axisymmetric body moves in K d , d > 2 along its symmetry 
axis, in a medium of chaotically moving particles; the medium is homogeneous, 
and distribution of the particles over velocities is the same at every point. The 
magnitude V of velocity is constant. The length h of the body along the axis is 
fixed, and the maximal cross section of the body by a hyperplane orthogonal to 
the axis is a unit (d — l)-dimcnsional ball. We consider the problem of finding 
the shape of body minimizing resistance of the medium. The main results of 
this paper are as follows. 

If d > 3, there are two different kinds of solutions. We shall describe them 
in the case d = 3; if d > 3, the description is quite similar. The solution of first 
kind is similar to the solution of classical Newton's problem, that is, its surface 
can be described in the same way as the surface of Newton's solution. The 
solution of second kind is a union of two bodies similar to Newton's solution 
"glued together" by rear parts of their surfaces. The length (along the direction 
of motion) of the front body is always more than the length of the rear body 
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(a) V = 1, h = 1.97 



(b) V = 1, h = 3.11 



Figure 1: Two solutions of the three-dimensional problem related to motion in 
a rarefied monatomic homogeneous ideal gas. The mean square velocity of gas 
molecules equals 1. 



turned over. The solution of first kind is realized for h < /i*, and of second 
kind, for h > h„, where h* = h*(V) > is a critical value depending on V. 
The function h*(V) goes to infinity as V — > +oo and to zero as V — > + . The 
examples of solutions of first and second kind are shown on Fig. 



and on 

Fig 1(b) respectively. Here and in what follows the body is supposed to move 
vertically upwards. 

If d = 2, the classification of solutions is somewhat more complicated. There 
are five different kinds of solutions: (a) a trapezium, (b) an isosceles triangle, (c) 
the union of a triangle and a trapezium, (d) the union of two is oscele s tria ngles, 
(e) the union of two triangles and a trapezium; see Fig. 2(a) - Fig 2(e) The 



solutions of first kind are realized for < h < vR, of second, for vP + < h < u*, 
of third, for < h < + u°_ , 
values 



iT!_ , of fourth and fifth, for h > 

u,(V), u°_ 



u 







These 



u°_(V) will be defined in section 4.2; 
one has < u° + (V) < ti*(V) < u*(V) + u°_(V). The solutions (a) - (d) are 
realized for any distribution of particles over velocities and for any positive V; 
the solution (e) is realized only for some special distributions and some values 
of V. The numerical computation of a solution of kind (e) is a hard task, which 
is unsolved as yet. 

In the limit cases, where the velocity of body is big and where it is small as 
compared to the mean velocity of particles, the shape of body of least resistance 
depends only on the length h, and does not depend on the distribution of par- 
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(c) h = 6 



(e) 
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(d) 
h 

7.83 



Figure 2: The two-dimensional problem. The solutions cor- 
responding to the cases (a) - (d) are calculated numerically, 
for the motion with velocity V = 1 in a gas; the gas param- 
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tides over velocities. In the first limit case the optimal body coincides with the 
solution of classical Newton's problem. In the second limit case, for d = 3, the 
optimal body is a second kind solution symmetric with respect to a plane per- 
pendicular to the symmetry axis, the inclination angle of the lateral surface at 
its upper and lower points with respect to this plane being 51.8°; and for d = 2, 
the optimal body is one of the four figures: (a) a trapezium if < h < 1.272; 
(b) an isosceles triangle if h = 1.272; (c) the union of an isosceles triangle and 
a trapezium if 1.272 < h < 2.544; (d) a rhombus if h > 2.544. In the cases (a) 
- (c) the inclination angle of lateral sides of these figures with respect to the 
base equals 51.8°, and in the case (d), exceeds this value. 

In a monatomic ideal gas the velocities of molecules are distributed according 
to Gaussian law. Suppose that the mean square velocity of molecules equals 1, 
then the kind of solution is determined by two parameters: velocity of the body 

V and its length h. We define numerically the regions on the parameter plane 
corresponding to different kinds of solutions; for some special values of param- 
eters we determine the shape of optimal body and calculate the corresponding 
resistance. This work is made in the two- and in the three-dimensional cases. 

Imagine an observer travelling with the body through the medium; he would 
detect a flux of particles falling on the immovable body. In fact, this picture is 
more convenient for us, and will be taken in the sequel. 

The paper is organized as follows. In section [21 the formulas for pressure of 
the flux on the body's surface and for resistance force are derived; two auxiliary 
lemmas of pressure distribution over the surface are formulated; and the prob- 
lem of minimal resistance is reduced to the form more adapted for studying. In 
sectional some auxiliary minimization problems are solved. Using these results, 
in the next section we solve the minimal resistance problem in general form. The 
solutions are different for the cases d = 2 and d > 3. In section [5J the obtained 
results are applied to the flux corresponding to a rarefied monatomic homoge- 
neous ideal gas of positive temperature. In appendix A, the auxiliary lemmas 
are proved, and in appendix B, asymptotic formulas for pressure functions as 

V — > + are obtained. 

2 Calculation of pressure and resistance 

2.1 Consider an immovable convex body B in Euclidean space M d , d > 2 and 
a flux of infinitesimal particles falling upon it. Velocities and masses of the 
particles in general are different. Let the function p(v) denote the distribution 
density over velocities of total mass of particles in a unit volume, so that for any 
two infinitesimal regions V, X C M. d having d-dimensional volumes |V| = dv, 
\X\ = dx, the total mass of particles that are contained in X and have velocities 
v G V equals p(v) dv dx. Therefore, the particles' distribution over velocities and 
density of the flux v — J Md p(v)dv are the same at each point and at each instant. 
It is supposed that v < oo. 

An individual particle hitting the body at x transmits the impulse m ■ 
2(v\n x )n x to the body, where m and v denote the particle's mass and its velocity 
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before the collision, n x means the outer unit normal to dB at x, and (• | •) means 
scalar product. Note that (v\n x ) < 0. 

Let S be an infinitesimal part of dB containing x, and let V be an in- 
finitesimal region containing v, of volume |V| = dv. The total mass of par- 
ticles, colliding with S in a time interval dt and having velocity v e V, equals 
dM = p(v) dv ■ (v\n x )~ \8\dt, where \8\ means (d— l)-dimensional area of 8, and 
z_ := max{-z, 0}. The total impulse transmitted by these particles equals 

dM ■ 2{v\n x )n x ~ —2p(v) dv ■ (v\n)_ 2 \S\dt ■ n x . 

Integrating this value with respect to v, one obtains the total impulse transmit- 
ted to S per time dt, 

{v\n x )J p(v) dv ■ \5\dt ■ n x . 

i. 

Dividing this value by \S\dt, one obtains that pressure of the flux at x equals 
Tr(n x ), where 

Tr(n) = -2 (v\n)_ 2 p(v) dv ■ n. (2.1) 



Integrating the pressure over dB, one gets the total force R(B) the flow is 
exerting on the body, 

R(B)= f n(n at )dH d - 1 (x), (2.2) 

JdB 



'dB 

jd-l 



where Ti means (d — l)-dimcnsional Hausdorff measure. 

2.2 Denote R+ := [0, +oo), and denote by Ad the set of functions a 6 
C 1 (1R+) such that the function a'(r)/r, r > is negative, bounded below, and 
monotone increasing, and 

/•run 

r 2 a{r)dr d < oo. (2.3) 



o 

Remark 1. Note that if <j\, U2 G Ad then a\ + ct 2 £ Ad- Also, if a, (3 > 0, 
a G Ad and d(r) = cta((3r), then a € Ad- 

From now on, we suppose that the density function p satisfies the condition 

A p(v) = a(\v + Ved\), where a £ Ad, V > 0, and is the dth 
coordinate vector. 

Note that the relation l|2.3|l implies that pressure ir(n) 12.1fl is always finite. 

Example 1. Consider a rarefied homogeneous monatomic ideal gas in M 3 of 
absolute temperature T > 0. The distribution density of molecules' mass over 
velocities equals o>j(|i;|), where 



3/2 mr 2 



G 



(the Maxwell distribution); here k is Boltzmann's constant and v is the gas 
density. Consider a body moving through the gas with constant velocity of 
magnitude V in the direction of third coordinate vector e^. In a frame of 
reference connected with the body the distribution density over velocities equals 
Ph(v) = 07i(|w + ^ e 3|)- It is easy to check that ah £ A3, so the condition A is 
fulfilled. 

Example 2. Let, now, a rarefied ideal gas of temperature T be a mixture of 
n homogeneous components, the ith component having density u% and being 
composed of monatomic molecules of mass rrii. Then the distribution density 
of molecules' mass over velocities equals cr„h,(|u|), where 

-A / mi \ 3 / 2 _mi£ 
i=l 

Taking account of remark 1, one concludes that a n h £ A^. As in the previous 
example, a body moves in the gas along the third coordinate axis with velocity 
V. In a frame of reference connected with the body the distribution density 
over velocities equals p n h{v) = <Jnh{\v + Ve^\) 1 therefore, the condition A is also 
satisfied. 

In the examples ^ and [21 the force of resistance of the gas to the body's 
motion is calculated according to l|2.1[l and l|2.2[) . 

In what follows, we shall also suppose the following condition to be fulfilled: 

B The body B is convex, compact, and symmetric with respect to 
the dth coordinate axis. Moreover, the maximal cross section of the 
body by a hyperplane orthogonal to the symmetry axis is a unit 
(d — l)-dimensional ball. 

By translation along the dth coordinate axis, the body can be reduced to 
the form 

B = {(x\x d ) : |x'| < 1, /_(|ar'|) < x d < -f+(\x'\)} , 

where a;' = (xi, . . . , x^—i), /+ and /_ are convex non-positive non-decreasing 
continuous functions defined on [0, 1]. The length h of body along the symmetry 
axis equals 

fc = -/+(0)-/_(0). 



Now, let us specify the formulas for pressure 7r(n) (|2.1(l and force R(B) H2.2(l . 
using the conditions A and B. 

At a regular point x + = (x\—f + (\x'\)) of the upper part of the boundary 
dB, the outer unit normal vector is 

"*+ = 1 ( /i(|z'|)^7, l) , (2.4) 

+ V/UM) 2 + i V +Kl u \x'\' J' 1 ' 
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and using (|2.1[) and taking into account axial symmetry of p with respect to the 
dth coordinate axis, one finds that pressure of the flux at this point equals 



7r(n x+ ) = ~P+ (/+0'|)) • n x+ , 



where 



p + (u) :-- 



1 



{u, 0,...,0, 1) 



(2.5) 



(2.6) 



. \/u 2 + 1 

Similarly, pressure of the flux at a regular point ir_ = (x\ /_(|a;'|)) of the lower 
part of dB equals 

where 



n(n x _) =v- (f'-(\x'\)) ■ n x _, 
1 



V/'(M) 2 + i v r - (M) i*"» 1 



and 



P-W 



1 



From l|2.6|l . (|2.9|l . and l|2.f|l one obtains 

(viU + £V d \ 2 



(u, 0, 



,o, -i; 



Pe(u) = £ 



1 + u 



p(v)dv, where e G { — , +}. 



(2.7) 
(2.8) 

(2.9) 
(2-10) 



Let us calculate R(B). The integral in the right hand side of l|2.2ll is the sum 
of two integrals corresponding to the upper and lower parts of dB. Changing 
the variable in each of these integrals and using the formulas (12 .4f) . 1)2. 5fl . 12. 7|) . 
and (12. 811. one obtains 



R(B) = / p+(J+(M)) • -/+(M)]3T, "I 

+ / P-(/i(k'l))- f/KM)]^, -i) 

i|x'|<i V F I / 

Next, using that the functions p £ (f^(\x , \)) are invariant, and the functions 
£ f'e(.\ x '\) iff are anti-invariant with respect to central symmetry x' — > —x\ one 
gets that 

hence 

R(B) 
Therefore 



Pe^dl'l)) • £ /^(|x'|) — T dl' = 0, EG {-, +}, 



|*'|<1 



p+(/;(| 3! '|))dx'+ / M/MM))^' e d . 

|x'|<l J|x'|<l / 



R(B) = -a d -! (Tl + (f + )+K-(f-))-e d , 
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where a^-i is the volume of a unit ball in M. d 1 , and 

pi 

■ d-l 



n e (f)= Pe (f(t))dt d -\ s g {-,+}. (2.ii) 



o 



Denote by M. (ft) the class of convex non-positive non-decreasing continuous 
functions / defined on [0, 1] such that /(0) = —ft. Note that any function / £ 
A4(ft) is differentiable everywhere except possibly on a finite or countable set, 
and /' is monotone, hence the integral (|2.11() is well defined for any / £ M.(h). 

2.3 Thus, the problem of minimal resistance takes the following form: 

minimize lZ+(f+) +72-_(/_) 
provided that / + and /_ are convex non-positive 
non-decreasing functions satisfying the relation 
-/+(0)-/_(0) = ft. 

It will be solved in two steps. First, given ft_ > 0, h + > 0, find 

inf Tl-(f) and inf TZ+(f). (2.12) 

feM(h.) yJ ' feM(h+) +VJ ' y ' 

Second, given solutions , f£ of the problems (|2.12|) . find 
R{k) := h+ $t= h ( n +(fh + ) + K-(fh_ 



2.4 Let us formulate two auxiliary lemmas. Their proofs are rather bulky, 
and are given in Appendix A. 

Lemma ^ states some properties of the functions p+ , p_ , which will be 
needed in the subsequent sections. 

Lemma 1. Let p satisfy the condition A. Then 

(a) there exist the limits Ym\ u ^ +00 p e {u) =:p e (+oo), besides 
p + {+<x>) +p_(+oo) = 0; 

(b) p E eC 1 (R+), andp' s (0) = ]im u ^ +oo p' e (u) = 0, e e {-,+}; 

(c) for u > 0, p' + (u) < p'_{u): 

(d) for u > 0, p' + (u) < 0, and for any u > 0, P-{u) > p_(+oo). 

Lemma 2 specifies the form of functions p + , p- for d = 2: the function p + 
has, in a sense, a simple behavior, and the behavior of p- may be complicated. 
This specification will be used in section 4.1 when constructing the body of least 
resistance in two dimensions. 

Designate a a 'P(r) = a(r) + a. a (fir), where a > 0, f3 > 0. By virtue of 
remark^] if a € ^2 then cj q,/3 S ^2, hence the function p a '@(v) — a a 'P(\v + 
Ve2\), V > 0, defined on R 2 , satisfies the condition A. Denote by p"' 13 , s £ 
{ — , +} the function corresponding to the density p a ^ , according to the formula 
H2.10f> . and denote by p"'P the maximal convex function defined on K + that 
does not exceed p"' 13 . 



9 



(a) 



(b) 



Figure 3: 



Lemma 2. Let d = 2. 

(a) If p satisfies the condition A 



then for some u > 0, p' + is monotone 



decreasing on [0, u] and monotone increasing on [u, +oo) (see Fig. 3(a)). 

(b) Suppose that a € A2, V > 0, and for any n > the function r n a(r) 
monotonically decreases, for r large enough. Then there exist a > 0, [3 > 
such that the set O a,/3 — {u : p°Z^ '(u) > p°_ (it)} has at least two connected 
components (see Fig. \3(b]\ O a ^ is shown bold-faced on the x-axis). 



3 Auxiliary minimization problems 

3.1 The following lemma reduces the minimization problems l|2.12ll to a simpler 
problem of minimization for a function depending on a parameter. 

Lemma 3. Let p e C(M+), d>2, A > 0, and let a function fh S A4(h) satisfy 
the condition 

(C\) fhO-) = 0j and for almost every t, u = f' h {t) is a solution of 
the problem 

t d - 2 p(u) + Xu -> min. (3.1) 
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Then f% is a solution of the minimization problem 

inf H(f), K(f)= [ p(f'(t))dt d - 1 . (3.2) 

Moreover, any other solution of fff.jjj) satisfies the condition (C\) with the same 
A. 

Proof. In fact, the problem 1)3. 2f> can be considered to be a degenerated case 
of the classical problem of optimal control [T2], an d the statement of lemma is 
a consequence of the Pontryagin maximum principle (see ^3]). The proof we 
give here, however, is quite elementary and does not appeal to the maximum 
principle (cf. |T5]V 

For any / S A4(h) one has 

t d - 2 P (f(t)) + Xf(t) > t d - 2 p{f' h {t)) + \f' h {t) (3.3) 

at almost every t. Integrating both sides of (|3.3|l over t E [0, 1], one gets 

^ f 1 P (f'(t)) dt*- 1 + A (/(I) - /(0)) > 

-d^rl! p{f ' h(t)) dtd ~ i + x Uh[i) hm ' (3 - 4) 

and using that /(l) < = /h(l) and /(0) = fh(0) = —h, one obtains that 
W) > K(h). 

Next, suppose that f £ M(h) and IZ(f) = IZ(fh), then, using the relation 
H3.4[) and the equality /(0) = //i(0), one gets that /(l) > fh(l) = 0, hence 
/(l) = 0. Therefore the inequality in l|3.4() becomes equality, which, in view of 
H3.3fl . implies that 

t d - 2 P {f'{t)) + Xf'(t) = t d - 2 P (f h (t)) + Xf h (t) 

for almost every t, hence u = f'(t) is also a solution of l]3.ip . on a set of full 
measure. Thus, / satisfies the condition (Ca). □ 

3.2 Assume, additionally, that p G C 1 (R+) and that p is bounded below. 
Denote by p the maximal convex function defined on K + that does not exceed 
p. The function p is also continuously diffcrentiablc, and for any h > 0, tt > 
one has 

p(u) > p(u) > p(h) + p\h) ■ (u - h). (3.5) 

Define the set O p :— {u : p{u) > p(u)}. Obviously, O p is open, and hence is a 
union of a finite or countable (maybe empty) set of disjoint open intervals. 

The following lemma specifies the solution fh of the minimization problem 
H3.2fl in the case d = 2. 
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Lemma 4. Suppose that the function p is bounded below, p € C 1 (R+), h > 0, 
< ; and denote by (h^\ /i^ + )) the maximal interval contained in O p such 
that h^' < h < ftW ^ ma?/ happen that = h = h^ + \ i.e., the interval is 
empty). 

Then fcW < +00, and i/ie following holds true: 
(a) The function fh defined by 

f h {t) = -h + ht, 

ifh^=h = h(+\ and by 

, f -h + ht-h if t<to, , , hW-h 

h{t) = \ , , ,(_). . . (+ w. .v ... ... where t = 



-h + h(-h + h^(t-t ) if t>t , 
is a solution of the minimization problem 

mj TZ(f), K(f)= [ p(f'(t))dt, (3.6) 

besides 

(b) inf K(f)=K(f h )=p(h). (3.7) 

fEM(h) 

(c) If f is a solution of J^3. 61) then at almost every t, the value u — f'(t) 
satisfies the relations p{u) = p{h) + p'(h) ■ (u — h) and u ^ O p . 

Proof, p is bounded below and p'(h) < 0, therefore there exists a value u > h 
such that p(u) = p(u), hence < u < +00. 

One has p(h^) = p(h^), p(h^) = p(/i (+) ), and 

p(/»C±)) = p{h) + p'{h) ■ (ftW - h). (3.8) 
From H3.5fl and l|3.8|l it follows that for any u 

p(u)-p(h^)>tf(h)-(u-h {±) ), 
and designating A = —p'{h), one obtains 

p{u) + Xu > p(/i (±) ) + Xh^K 

This means that both hS~) and h^ minimize the function p(u) + Xu. 

Further, one easily sees that fh 6 M.{h), = 0, and the function f' h 

takes the values h^~) and (which may coincide). Applying lemma 01 one 
obtains that fh is a solution of the problem (|3.6|l . 

If ft(-) ^ then 

ft(/fc) = J° p(h ( --1)dt + J p{h^)dt = 

° hl+ '-" p(^) + JAp(^- 



h(+) - hi-) ^ y ' - ft(-) 
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On the other hand, excluding p'(h) from the relation (|3.8|l . one obtains 



hW-hi-) ry ' hW-hi-) 

hence K(f h ) = p(h). If = = h then TZ(f h ) = p(h) = p{h), so the 
formula l|3.7|) is true. 

Let, now, / be a solution of l|3.6[) . By lemma for almost every t the value 
u = f'(i) minimizes the function p(u) + Ait, hence p(u) + Xu = p{h^) + Xh.( + \ 
and substituting A = —p'{h) and using that p(ft/ + ') = p{h) + p'(h) ■ (h^ — h), 
one obtains 

p(u) =p(h) +p'(h) ■ (u-h). 
Taking into account 13. 5L one gets that 

p(u) =p(u) =p(h) +p'(h) ■ {u- h), 

hence u ^ O p . Lemma 0] is proved. □ 

3.3 Suppose, in addition to the previous assumptions, that there exists the 
limit lim u ^+oo p(u) =: p(+oo) and that for any u G R+, p(u) > p(+oo). Then 
p'(u) < for any u. Denote B — — p'(0); one has B > 0; the function p' is 
continuous and monotone non-decreasing from B to 0. 

Suppose that d > 3, and introduce an auxiliary notation: lu = -335, 
q(u) = Q(u) = f q(v) dv. Both function q and Q are continu- 

ous and monotone non-decreasing on R + ; q changes from B~^ to +oo, and Q 
changes from to +oo. 

Lemma 5. Let d > 3. h > 0, p G C 1 (IR + ), and let th ere exist the limit 
lim l[ _ ++00 p(u) — p(+oo) and for any u G R + , p(u) > p(+oo). Then 

(a) the set of values U > satisfying the equation 

is non-empty. 

(b) Let U be a solution of 13. ,9|) and to = -^jj- Introduce the function fh(t), 
t G [0, 1] as follows: 

forte[0,to], f h (t) = -h; 

for t G [to , 1] , fh is defined parametrically, 

uq(u) - Q(u) 
h = -h+ W) , (3.10) 

t=^j- y ue[0,U}. (3.11) 
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The function fh is defined correctly, is strictly convex on [to, 1], and is a unique 
solution of the minimization problem \3.fy) . 
(c) The minimal value of 1Z equals 

inf K(f) = K{f h ) = p(U) + -7)7(^7 ■ (3.12) 

Proof, (a) For arbitrary c > 0, one has 

q(U) J \ q(U)J -Jo \ q(U)J 

and taking into account that lim;y_^ +00 q(U) = +00, one concludes that for U 
sufficiently large, U — > c/2. This implies that the continuous function 

U — goes to +00 as U — > +00; besides it vanishes at U = 0, hence the set 
of solutions of l|3.9[l is non-empty. Therefore, (a) is proved. 

Denote by S(h) the set of points u such that p{u) = p(h) + p'(h) ■ (u — h). 
Note that S(h) coincides with the connected component of O p containing h, if 
h € O p , and S(h) = {h} otherwise. Thus, S(h) is a closed segment containing h; 
two segments S(hi), S(h2) either coincide or are disjoint. The condition p(u) > 
p(+co) implies that all segments S(u) are bounded. Obviously, the family of 
non-degenerated segments (i.e., of those that are not singletons) is at most 
countable. Denote by S the union of non-degenerated segments. The function 
p'(u) is monotone increasing on R + \ S and is constant on any non-degenerated 
segment S(u) C S, hence the set {p'(u), u £ S} is at most countable. 

Let < ui < it 2 . After simple algebra one obtains 

Q(ui)~ 
1 — 1 — ^ ~ 



Q(ui) 



q(ui) q(u 2 ) 



(q(u 2 ) - q{v))dv. 



(3.13) 



Both terms in the right hand side of (|3.13|l are non-negative, hence the function 
U — ^an is monotone non-decreasing. If both ui and u 2 are solutions of (|3.9|l 
then both terms in H3.13|l are equal to zero, which implies that q is constant on 
[til, U2]; that is, p' is constant on [u±, u 2 ]', or, equivalently, u\ G S(u 2 ). This 
implies that the solution set of (|3.9() coincides with a segment S(u). 

(b) The relations H3.10|) and l|3.11|l define continuous functions fh and t of u, 
varying from — h to and from to to 1, respectively, when u passes the interval 
[0, U\. Moreover, the function t = q(u)/q(U) is monotone non-decreasing, each 
set {u : q(u)/q(U) — t} coincides with some segment S(u), and fh is constant 
on any such segment. This means that the function fh{t) is defined correctly. 

Let us calculate the left-hand and right-hand derivatives /^(i - ), 



Designate by [u (t), u + (t)] the interval {u : q{u)/q{U) — t} and put u 



'(*). 
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Let the values fh + Afh, t + At correspond to the argument u + Au, and At > 0. 
One has 

q(u + An) — q{u) 



Afi, 



M - HP) • 

(u + Au) q(u + Au) — Q(u + Au) u q(u) — Q(u) 

: W) W) 

(q(u + Au) - q(u)) + C +Au (q(u + Au) - q{v))dv 



llGIlCC 

M = „ + r A "f + f>-«;y (3.i4) 

At J u q(u + Au)-q(u) 

The integrand in the right hand side of i|3.14[) is less than 1; due to definition 
of u = u + (t), one has Au — * + as At — > + , therefore 

&(t + )= lim M = „+ (i ). 

At->0+ At 

Similarly, one finds 

/i(t") =«"(*). 

Both functions u~(t) and u + (t) are positive, and for any t% < t^ one has 
u~(t\) < u + (t\) < u~(t2) < u + (t2). Therefore, the function fh is mono- 
tone increasing and strictly convex on [to, 1]; moreover, it is constant on [0, to], 
f h (0) = -h, and /(l) = 0. Thus, it is proved that f h G M(h). 

For any t G [to, 1], except possibly a countable set of values, one has u~(t) — 
u + (t) '.— u G R+ \ S, hence there exists the derivative f' h (t) = u. For any u ^ u 
one has 

p(u) > p{u) + p'(u) ■ (u — ll), 
and using that p(u) > p(it), = one obtains 

p(u) > p{u) + p'{u) ■ (u — u), 

hence 

p(u) — p'(u) ■ u > p(u) — p'(u) ■ u. (3.15) 

Recall that t = = and w = tV One has t fi - 2 = ffi, and 

multiplying both parts of (|3.15|) by t d ~ 2 and designating —p'(U) — A, one 
obtains that 

t d ~ 2 p(u) + \u> t d ~ 2 p(u) + Xu 

for any u ^= u. Thus, u = f' h {t) is a unique value minimizing the function 
t d - 2 p(u) + Xu. 

Let, now, t G (0, to). For u > one has 

p(u)>p(0)+p'(0)-u. 
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Using that p(u) > p(u), p(0) = p(0), ^ = p'(0)/p'(U) = -p'(0)/A, one 
obtains 

P(«) > P(0) +p'(0) • « = p(0) - At^u, 
hence for any t £ (0, to) 

p(u) + \t 2 - d u > p(0), 

therefore the value flit) = is a unique minimum of the function t 2 p(u)+\u. 
Applying lemma one concludes that fh is a unique solution of Q3.2j l. 
(c) One has 

Tl(fh) = f° Pi^df 1 - 1 + [ p{f h {t))dt d -\ (3.16) 



The first integral in the right hand side of (|3.16|) equals 



1{U) 

Denote U = infS(U). Changing the variable in the second integral t = q(u)/q(U), 
u € [0, U] and taking into account that for almost every t, f' h (t) = u, one ob- 
tains that the second integral equals 



g( M ) 



dp{u). 



Summing the first and the second integrals and taking into account that q{U) = 
q{U) and d— 1 = 1 + l/u), one obtains 

1Z(f h )=p(U) - [ U " dp(u). (3.17) 



The integral in 1)3. 17|) can be represented as the sum 
r U 

(•■■)= 

l[0,U]\S T JS - 



JO J\0.U]\S , JSi 



where Si = S(ui) are non-degenerated segments whose union gives [0, U] H S. 
If u e [0, U] \ S, one has p'(u) = p'(u) = — q^)^ 1 ^ , hence 

f / \ / , 

i[o,&]\5 l "' J y[o,&]\ S g(t/) 1+1 /- 

Next, taking into account that the function q is constant on Si and that at 
cndpoints of Si, p and p coincide, one obtains 

{ - )= JsXW)) dU= -Js,Wf^ dU - 



Si 
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Summing these integrals, one gets 



U / I \ \ 1+1/ V f U I \ 

qiUA Mu) = - ^MlZ^ (3.18) 



MU)J ^ ' Jo q(u) 1+1/ 

Further, the function q(u) = is constant on [U, U], hence 



g( J )1+1 /. du = -i(ur 1/uJ (u -u) = p'(u) (u-u)= p(u) - p(u). 

( 3 - 19 ) 

Using that p(U) = p(U) and applying lj3T7f) . (|3~T%|) . and (|XT^|i . one gets 

and recalling that Q is the primitive of q and 1 + l/u> = d — 1, one comes to the 
formula l{3TiSj) . □ 

4 Solution of the minimal resistance problem 

4.1 Two-dimensional problem 
4.1.1 Minimization of 1Z+ 

From statement (a) of lemma |21 it follows that there exist values u\ > and 
B+ > such that 

P+ «)- P+ (0) 
o = = -B+ 



P+( u ) 



and 

p+(0) ~B + u if < u < u° + , 
p+(u) if u > m?.. 

This implies that O p+ = (0, u ^). Applying lemma one obtains that there 
exists a unique solution f£ of the minimization problem 



/EM(d) 

defined by the relations 



inf ft+Cf), K + (f)=[ 1 p + (f 
M(h) J 



(t))dt, 



f+u\ J ~ h forte [0, t ], a ,s 

j \ + < -(i-io) forte[t ,l], 1 j 

to = 
if < h < u+, and 

/+(*) = -fc + ht, 
if h > ir}_. The minimal resistance equals 

mf ■ n + (f) = n + (f+)=p + (h). 

feM(h) 
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4.1.2 Minimization of 11- 

Note that p'_(0) = and p'_(0) < 0, hence O p _ contains an interval (0, u°_), 
u°_ > 0, besides P-(u°_) = p_(0) + pL(0) ■ u°_. Denote 5_ = -p'_(0) and 
represent the open set O p _ as the union of its connected components d = 
(u~ , uf), O p _ = UjOj. We shall suppose that the set of indices {i} contains 1 
and that 0\ — (0, u°_ ). Statement (b) of lemma[21and example [2] imply that in 
some cases (for example, when considering pressure distribution of a mixture of 
two homogeneous rarefied gases on the rear part of surface of a moving body) 
O p _ has at least two connected components. 
Consider the minimization problem 



inf 

f<SM(h) 



K-(J), K-(J)= [ P- (/'(*))*• (4-2) 
Jo 



Applying lemma ^ one obtains that there exists a solution f h of this problem, 
besides 



For < h < u°_ one has 

= { -h + u°_ -(t-to), if *> *o, (4 ' 3) 
t = l-h/u°_. 

For fteR\ O p „ one has 

fh{t) = -h + ht. 
Finally, for h £ Oi, i ^ 1 one has 

f -h + u^t, if t<ti , , 

7/1 W 1 -fc + u^ + «+(*-*,), if * > *< , 1 j 



/; 



Notice that //j needs not be the unique solution of 1)4. 2|l . In some degenerated 
cases it may happen that the right endpoint of some interval Oi coincides with 
the left endpoint of another interval, uf = uj , i ^ j\ then there exists a 
continuous family of solutions of l|4.2[) : the derivative of any function from this 
family takes the values uj , uf, and Uj. 

4.1.3 Solution of the two-dimensional problem 

Thus, the problem of finding 



R(h) = h jnf_ =h (ll + (f+ + ) ) 
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amounts to the problem 



min ph(z), where p h (z) = p + (z) + p-(h - z). (4.5) 

0<z<h 

The functions p'_(u), p' + {u) are continuous and monotone non-decreasing, 
hence the function p' h (z), < z < h is also continuous and monotone non- 
decreasing. 

Using statement (c) of lemma ^ one concludes that B + > Indeed, if 
u°_ < u+ then 

p_(u°_) - p_(0) p+(u°_) - p+(0) p+(u°_) - P+ (Q) 

-E>- o > o o = 

u_ u_ u_ 

and if u°_ > then 

-B- = pUu -) > p' + {u°_) > MO = -B+- 

Thus, there exists a unique value it* > vP + such that p', (it*) = p' + (u*) = —B-. 
Consider four cases: 

1) < h< u° + ; 

2) u+ < h < u*; 

3) < h < it* + 

4) h > + u°_ . 

In the cases 1) and 2), for < z < h, one has p' h (z) < + B- = 0, 

hence z = h is a unique value of argument minimizing ph- Therefore, the optimal 
values of h + and h— are h + = h, h— = 0, and fZ _ = 0. 

1) < h < The function f^ +=h is given by l|4.1|) . The body of least 
resistance is a trapezium, the tangent of slope of its lateral sides being equal to 



(see Fig. 2(a) i. The minimal resistance equals 



R(h) = K+ (/+ ) + K- (4" ) = P+ (0) - B + h + p_ (0). 
2) w+ < h < u*. Here one has f£ = ) l {t) = —h + ht, hence the optimal body 



is an isosceles triangle (see Fig. 2(b) I , and 



R(h)=p+(h)+p-(0). 

In the cases 3) and 4) one has p' + (h) > —B-, hence p' h (h) = p' + (h)—p'_(0) > 
0. On the other hand, p' h {u a + ) = p' + {u a + ) - p'_(h - u\) < -B + + B_ < 0. 
Moreover, using statement (a) of lemma|3 one finds that the function p' + (u) = 
p',(u), u £ [u® + , h] is monotone increasing, hence p' h is also monotone increasing 
on this interval. It follows that the function p^ has a unique minimum z € 
«, h) and /+ =2 (f) =-z + zt. 

3) < h < One has p' h { u *) = p' + (u*)~p'_(h—u*) = —B-+B- = 0, 

therefore ph reaches its minimal value at z — u*; thus, the optimal values of h+ 
and h- arc h + = u*, h- = h — u*. The function fr =h _ u is given by 14.3(1 . 
Here the optimal body is the union of a triangle and a trapezium, as shown on 
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Fig. 2(c) The tangent of slope of lateral sides of the trapezium equals 



The minimal resistance equals 

R(h) = p+{u*) +p-{0) -B_(h- u*). 

4) h > + u°_ . One has p' h (h - u°_) = p' + (h - u°_) + B_ > 0, hence the 
minimum of ph is reached at a point z € (ir], ft. — it?.], and the optimal values 

= z, h— = h — z, as well as the minimal resistance, are obtained from the 
relations 

h + + h- = h, 
p' + (h + )=p>_(h-), 
h+ > u\, /i_ > u°_, 
R(h)=p+(h+)+p-(h-). 

Here one should distinguish between two cases. 

4a) If /i_ 6l \ O p _ then fr (t) = — h- + h-t, and the optimal body is a 
union of two isosceles triangles with common base, of heights h + and h- (see 



Fig. 2(d) I 



4b) If belongs to some interval Oi = (u. , u^), i ^ 1, then /. is given 
by 1)4.4(1 . and the optimal body is the union of two isosceles triangles and a 



trapezium (see Fig. 2(c) I 



Note that the case 4b) is realized for the values h from an open (maybe 
empty) set contained in (u* + u°_, +oo). This set is defined by the parameters 
a and V. The case 4a) is realized for the values h from the complement of this 
set in (ii* + u_, +oo), which is always non-empty. 

4.2 The problem in three and more dimensions 

Let d > 3. Using lemma0 one concludes that there exists a unique solution /| 
of the problem 

inf TZ e (f), TZ e (f) = f Pe (f(t))dt d -\ 

besides 

q E (U) d 1 

where U is defined (not necessarily uniquely) by the relation 

Qe{U) 

here q e (U) = W^UT 1 /^, Q £ (U) = tfq s (u)du. 
Thus, the problem 
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amounts to the following problem: 



where 



n+ (ii_|_)+n_ (u_ J — /i 



/ \ - / x . Q+(«) , , _ / \ Q-(«) 

-W=P+(«) + — ? w i ) r_(u)=p_(u) + 



and 



_v ; r-\ j < q _(u)d-i 

h+(u) = u T-^r, h-(u) = u -— , u > (J. 

g + (u) g_(tt) 

The functions r e and p e , £ G {— , +} are monotone non-increasing, and h e is 
monotone non-decreasing from to +oo when u € R+, besides any interval of 
constancy of one of these functions is at the same time the interval of constancy 
of two others. For each z > choose u such that h e (u) = z and put t^ e '{z) := 
r s (u), n^(z) := p' £ (u). From the stated above it follows that the functions r^ £ ) 
and 7r( £ ) well defined on R + and are monotone decreasing. Denote 

r h ( z ) = r (+ )(z) +^-\h-z). 

After some algebra one obtains that the function th is differentiable and 

r' h (z) = (d- 1) (p' + (u + )-p'_(u-)), (4.6) 

where the values ti+, U- are chosen from the relations h + (w + ) = z, h_(it_) = 
h — z. Both values in the right hand side of i|4.6[l . p' + {u + ) = 7r( + )(z) and 
— = — 7P — z), are monotone increasing functions of z, hence r^(z) 
is also monotone increasing from r' h (0) = (<f — 1) (p+(0) — p'_(U-)) to rj^/i) = 
(<i — 1) (p' + (U+) — p'_(0)), where U+ and U- are defined from the relations 
h+(U+) = h, h_(Z7_) = h. Note that p' + (0) = -B+ and p'_(U-) > -B_ > 
—B + , therefore r^(0) < 0. 

Recall that it* is defined in the subsection 4.1.3 by p' + (u^,) = — £>_. Designate 

K :=h+(tt*) = - B 1 ^ Q + {u*) (4.7) 

and consider two cases. 

1) h < /i*. One has h+([/+) < h+(u*), hence E/+ < u*, therefore r' h (/i) = 
(d-l)(p' + (U + ) + B_) < (d-l){p + {u*)+BJ) = 0. This implies that r^(z) <0 
for z G [0, h), hence the function r^ has a unique minimum at the point z = h, 
which corresponds to the values h + = h, h— = 0. The minimal resistance equals 

R(h) =p+(u + ) + Q + {u+)q + {u + )- d+1 +p-(0). 

2) h > K. One has U+ > u„, therefore r' h (ti) = (d- 1) (£/+) + B_) > 0. 
On the other hand, rl(0) < 0. Hence, there exists a unique value z S (0, h) such 
that r^(z) = 0. Thus, the function r^ has a unique minimum at z; the optimal 
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values of h + , h- arc h+ = z > 0, h- = h — z > 0. These values and the related 
auxiliary values U- , u + are uniquely defined from the system of four equations 

h+ = u+ - Q + (u+)/q + (u+), 
h- = u- — Q-(u-)/q-(u-), 
h + + h- = h, 
P'+{u+) =p'_{u-), 

and the minimal resistance equals 

R(h) = p + (u + ) + Q + (u + )q + (u + r d+1 + p_(«-) + g_(w_) g _( U _)- d+1 . 

4.3 The limit cases 

Consider heuristically the limit behavior of solutions as V — > +oo and as V — > 0, 
with fixed /i and a. We shall denote the pressure and resistance functions by 
p±(u,V) and R(ft, V), thus explicitly indicating dependence of these functions 
on the parameter V. 

4.3.1 V -» +oo 

Denote by p e (u, V) = V~ 2 p e (u, V), e G { — , +} the reduced pressure, and by 
R(h, V) = V~ 2 R(h, V), the minimal reduced resistance. One has 

p+(u,V) = l/(l + u 2 ) +o(l), 
P-(u,V)=o(l), p'_{u,V) = o{l), V^+™ ; 

in other words, as V — > +oo, the functions V") and V) tend to 

1/(1 + u 2 ) and to 0, respectively. These limit functions determine pressure 
distribution on the front part and on the rear part of body's surface in Newton's 
classical problem. 

Consider the cases d — 2 and d — 3 separately. 

d = 2 If h < 1 then for V sufficiently large, the figure of least resistance 
is a trapezium, and the inclination angle of its lateral sides tends to 45° as 
V — ► +oo. If h > 1 then for V sufficiently large, the figure of least resistance is 
an isosceles triangle coincident with the solution of two-dimensional Newton's 
problem. 

d = 3 For V sufficiently large, the body of least resistance is the first kind 
solution. The front part of its surface is the the graph of a function defined on a 
unit circle; as V — > +oo, this function uniformly converges to the function that 
describes Newton's classical solution with the same h. 

The case d > 3 is similar to the three-dimensional one. 

Finally, the limit value of minimal reduced resistance R(/i, oo) is equal to the 
resistance of Newton's optimal solution multiplied by the density of particles' 
flow v = lMicr{\v\)dv. 
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4.3.2 F->0+ 

In this limit case one has 

p E (u,V) = ebW + V , +o(n {-,+}, 

V 1 + vr 

p+OO r+OO 

6 (2) = - / a(r)r 3 dr, c (2) = 4 / <r(r)r 2 dr, (4.8) 
2 Jo Jo 

6 (3) = — CT ( r ) r 4 ^ c (3) = 2vr / cr(r)r 3 dr. (4.9) 

3 Jo Jo 

This formula will be derived in Appendix B. One readily obtains that vP_ 



where 



and u° + tend to the value a := J (1 + \/5)/2 » 1.272, and B± = V- a" 5 + o(V). 
Taking into account that p' + (u) < p'_(u) < 0, one concludes that tends to 
the same value a, and u+ + u* tends to 2a. 

Let us describe the shape of optimal body and determine the minimal re- 
duced resistance K(h, V) — y _1 R(/i, V) in the limit V — > + . We shall distin- 
guish between two cases: d — 2 and d = 3. 

d = 2 (a) < /i < a: the optimal body is a trapezium. 

(b) h = a: an isosceles triangle. 

(c) a < h < 2a: the union of a triangle and a trapezium. 

(d) h > 2a: a rhombus. 

In the first three cases, the tangent of slope of lateral sides of optimal figures 
equals arctana w 51.8°, and in the last case, exceeds this value. The examples 
of optimal figures are shown on Fig. 

The minimal reduced resistance equals 



where 



R(h, V) = 2c {2) p(h/2) + o(l), V -> 0+, (4.10) 



, 1 - a 5 u, if u < a, 

1 1 /VI if M>0. 



d = 3 Let 



Q(u) 



if u < a 
if u> a, 



a 5 u if u < a 



and let U be a (unique) solution of H3.9(l . Define the function like in the 
statement (b) of lemma |SJ Then the body of least resistance is 

{(x\x 3 ) £l 3 : M<1, \x 3 \ < -f h (\x'\)} , 

where x' = (xi, x^). Thus, the body is symmetric with respect to the horizontal 
plane {23 = 0}. The front and rear parts of its surface contain flat circular disks 
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(a) h = 0.5 



(b) h = 1.27 




(c) h = 1.77 (d) h = 3 



Figure 4: Two-dimensional case. Solutions in the limit V — > + are shown 
by solid line. The corresponding solutions of Newton's classical problem are 
shown by dashed line. The case (b) is the unique one where these two solutions 
coincide. 
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Figure 5: Three-dimensional case, h = 1. The solution in the limit V — > + 
(solid line) and the solution of Newton's classical problem (dashed line) . 

of equal radius, and the angle of slope of lateral surface near these disks equals 
arctana « 51.8°. 

On Fig. 03 the projections of two optimal bodies of height h = 1 on the 
plane Ox\x^ are shown. The bodies are: the body of least resistance in the 
limit V — > + and the solution of Newton's classical problem. 

The minimal reduced resistance equals 

where p is defined in (|4 . 1 II) , and U is given by (|3.9I) . 

It is interesting to note that in these limit cases, the shape of optimal body 
does not depend on the distribution a; moreover, the reduced minimal resistance 
is proportional to v in the limit V — > +oo, and is proportional to the factor cW 
given by l|4.8|l . (|4.9|l in the limit V — ► + . This factor can be interpreted as the 
sum of absolute values of impulses of particles of the medium in unit volume, 
in the frame of reference associated with the medium. 

5 Gaussian distribution of velocities: exact so- 
lutions 

Suppose that the function p = py is the density of circular gaussian distribution 
with mean —Ved and variance 1, i.e., 

pv (v) = a(\v + Ve d \), where a{r) = (2tt)-^ 2 e^' 2 . (5.1) 

This function describes the particles' distribution over velocities in a frame of 
reference moving in a homogeneous monatomic ideal gas, where the velocity 
of motion equals V and the mean square velocity of molecules equals 1 (see 
example . 

The function py satisfies the condition A, therefore the results obtained in 
the previous section can be applied in this case. Here, the pressure functions 
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p±{u,V) are calculated analytically in the cases d = 2 and d = 3, and then, 
using numerical simulation, the following results are obtained: 

1) The parameter set V-h is divided into several subsets corresponding to 
different kinds of solutions. This partition is shown on Figures^! and [5] 

2) The least resistance K(h, V) is calculated for various values of h, V. The 
results are shown on FiguresT71 and ITU1 

3) For several values of parameters h and V, the body of least resistance is 
constructed. Two such bodies are shown on Figures^] 0(a)-(d) and|H] 

Here the value V is allowed to vary, so the pressure and resistance functions 
are designated by p±(u, V) and R(h, V) instead of p±(u) and R(h). 
Consider the cases d — 2 and d = 3. 



5.1 Two-dimensional case 

Fixing the sign "+" and passing to polar coordinates v = (— r sin 95, —r cos <p) in 
the formula 1)2.10(1 . one obtains 

f f r 2 (cos w + u sin w) , 2 . , 
p+(u, V) = J J -i \ + u2 ^+ p+(r, <p, V) rdrdip, (5.2) 

where z+ := max{0, z}, and p+(r,ip,V) is the gaussian density function pv 
1)5. ljl written in the introduced polar coordinates, 

p + (r > <p,V) = ±e-^ 2 - 2Vr ™*+ v2 \ (5.3) 

Next, fixing the sign "— " and introducing polar coordinates in a slightly different 
manner, v = (— r sin ip,r cos<p), one obtains 

f f f 2 (cos w + u sin w) , 2 , 
p_ (u,V) = -jj -± ^ + ^ 2 * >+ p- (r, V) rdrdp, (5.4) 

Here p_ (r, ip, V) is the same density function pv 15. 1|) written in these coordi- 
nates, 

P-(r,<p,V) = _L e -^(- 2 +2^cos^)^ (5 5) 

27T 

Combining the formulas 1)5.2)1 . 15. 3)1 . 1)5.4)1 . and 1)5. 5|l . one comes to the more 
general expression 

/ e~ y2 / 2 /" /" (cosgj + usin</?) 2 _i r 2 , 9 v , 
p e (u,l/)=e^— // ^ \ + u 2 e r drd ^ ' 

cos (p+u sin (/>>0 

where e G {— , +}. Passing to the iterated integral and integrating over r, one 
obtains 

e~ v I 2 f (cos(fi + usina?) 2 , . , 

p E (u,V)=e / i V 2 l(eV cos ip) dtp, 5.6 

cos i^+u sin (/?>0 
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Figure 6: Two-dimensional case. Four regions shown on the parameter space 
correspond to four kinds of solutions. 



where 



and erf(x) = ^= J^e * eft. Changing the variable r = ip — arcsin(u/\/l + tt 2 ), 
one finally comes to 



Pe{u,V) = £ 



e y2/2 r /2 , , / _ cos t — u sin t . , 
cos^ t I \ eV , — dr. 



tt/2 



vT 



(5.7) 



Numerical simulations are done using Maple and verified by Matlab. 
Graphs of the functions 

h = u° + (V), h = u*(V), h = u*(V) + u°_(V), 

are shown on Figure^ where the values vR_, uP_ and it* (which are the functions 
of V) are defined in subsections 4.1.1, 4.1.2 and 4.1.3, respectively. These graphs 
separate the parameter space into four regions corresponding to the four 
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different kinds of solutions. The lower function h = it+(V) tends to 1 as V — ► oo. 
Further, at V — 0, the lower, the middle, and the upper functions take the values 

a, a, and 2a, respectively, where a = y (1 + V^)/2 w 1.272: limy^ M +(^0 = 
a, limv^o u *(^) = a , limy^o(w*(V) + u°_(V)) = 2a. 

The solutions of fifth kind (union of two triangles and a trapezium) were not 
found in numerical simulations. (These solutions correspond to the case where 
the set O p _ t v = {u : p-(u, V) < p-(u,V)} contains at least two connected 
components.) We believe that in the considered case corresponding to the gaus- 
sian distribution py this kind of solutions does not appear at all. (Notice that, 
according to the statement (b) of lemma this kind of solutions does appear 
for some distributions corresponding to mixtures of homogeneous gases.) 

Further, using the formulas from subsections 4.1.1 and 4.1.2, the func- 
tions / + and /_ are calculated, which allow one to construct the optimal fig- 
ures (Fig.|21(a)-(d)), and using the formulas from 4.1.3, the minimal resistance 
R(h, V) is calculated. The graphs of reduced minimal resistance R(/i, V) — 



V 2 R(/i, V) versus h are shown on figures 7(a) and 7(b) for several values of V . 



5.2 Three-dimensional case 

Fix the sign "+". In spherical coordinates v — (— r ship cos 9, — rsin<^sin#, —rcosp), 
r > 0, < <p < 7r, — 7r < 9 < 7r the formula H2.10|) takes the form 

f f f r 2 (cosp + usmpcos9), 2 , N „ , , , A 

p+(u,V) = JJJ -i ^ X + J ^ p+(r,p,9,V)r 2 smpdrdpd9, 

where 

p + (r,p,9,V) = -^±^e-i( r2 - 2V ^+V 2 ). 

Now, fix the sign "— ". In spherical coordinates v = (—r smp cos9, — rshxpshx6, rcosp), 
one has 

f f f r 2 (cos p + u sin p cos 9), 2 , , 
P-(u,V) = -JJj g 1 + u * ^ p- (r, p, 9, V) r 2 sin p drdpd9, 

where 



1 

(2tt) 3 / 2 

Summarizing, one comes to the formula 



p e (u,V) = e 



p_(r, p, 9, V) = e -|(r= + 2^cos V+ ^)_ 

es to the formula 

e~ v I 2 f f f ( cos p + u sin </? cos&Y 



(2tt) 3 / 2 77 7 1 + w 2 

cos v?+it sin tp cos #>0 



. e -|r 2 + £ Vrco S¥3r 4 sin ^ rfr .^ 6/) £ e {_ +}. ( 5 . 8 ) 

We shall use two formulas, which are easy to verify. First, 



1„2 



e 2 



r 2 +eVr cos ^4 ^ = j^y CQf , ^ ^ 



2N 
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where 



I{z) = e z2/2 (3 + 6z 2 + z 4 )(l + ed(z/y/2)) + 5z + z 3 



Second, 



where 



and 



(cos ip + u sin y> cos 9) z 



cos sin cos #>0 



1 + u 2 



if 



d0 



J (it, cost/?), 



1 < C < -u/Vl + ' 
J(u,C) = { Ji(u,0, if ICI <u/VT + l? 

J 2 (m,C), if u/VT+u* < C < 1, 



j 1 (u,0 = T -r-2 (2C 2 + « 2 (i - C 2 )) + 3CVV - C 2 (i + « 2 ) 

1 + u z L 

c 



l + u 2 



[2^+^(1-0], 



?o = arccos — - 



Taking into account these formulas and changing the variable £ = cos tp in l|5.8|) , 
one gets 



p e (u,V) = e 



(2tt) 3 / 2 



I(eVt)J(u,() dC 



e- v / 2 
(2tt) 3 / 2 



tl/Vl + M 2 



7(^0 -M^C) 



J(eVC)^(«,0dC 



Next, one numerically calculates the function h*(V) according to the formula 
(|4.7|l . This function is shown on Fig. 9(a) it divides the parameter set Ml into 
two subsets corresponding to two different kinds of solutions. The function 
/i*(V) looks like linear, but is not; the graph of its derivative is shown on 
Fig. 9(b) The function K(h, V) is calculated according to the formulas given in 
subsection 4.2; the graphs of K(h, V) versus h are shown on Fig. ^| for several 
values of V. On Fig.[5]the examples of solutions of the first and the second kind 
are presented, for parameters indicated there. 



Appendix A 
Proof of lemma [T] 

Changing the variable v — rv, r > 0, v € in the integral H2.10|) and using 

condition A, one obtains 



Pe(u) = e 



( dn d -\ v ) ( 

Js- 1 - 1 Jq 



2 <>1 M + £v d)J , > d _! , 

r ; — : — ^ P\ rv ) r dr 



l + u 2 
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(a) V = 1, h = 1.97 (b) V = 1, h = 3.11 

Figure 8: The solutions in three-dimensional case, for the distribution function 
pv H5.ip . 



31 



0.5 1 1.5 2 2.5 3 3.5 

V 



0.5 1 1.5 2 2.5 3 3.5 

V 



(a) h,(V) 



(b) dh t /dV 



Figure 9: Three-dimensional case. The function h*(V) divides the parameter 
plane V-h into two subsets corresponding to the two kinds of solutions. 
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(yru + ev d ) 2 
l + u 2 



S d-i 



where 

,.+00 (-+OC 

p{v) := I r 2 p{rv)r d ~ l dr = / r 2 a(y/r 2 + 2rV v d + V 2 ) v d ~ x dr. 
Jo Jo 

Substituting u = tan ip, ip g [0, 7r/2], one obtains 

p e (tanip) — e I {vi sinip + ev d cos<p)_ 2 p(v) dH d ^ 1 {v). 



Substitute "+" for e and consider the rotation T v that sends the vector 
(sin ip, 0, ... ,0, cos <p) to e d and leaves the vectors , i = 2, . . . , d— 1 unchanged. 
For any v G M d one has T v v — (cos<pv\ — sin <pv d , ^2, ■ • ■ , y d-\i sinyz'i + 
cosyi^)- Changing the variable T v v = lo, one gets 

p+(taa<p)= f LujpiT-'^dH*- 1 ^), (A.l) 
■/si -1 

where S^ _1 := {lu e S"*- 1 : w d < 0}. Designate 

g(z):= r 2 a(^r 2 + 2rVz + V 2 ) r d ~ x dr, \z\<l; (A.2) 
Jo 

obviously, p(y) = g{vd)- Using condition A, one concludes that the function g 
is continuously differentiable, and its derivative 



, f°° 2 <r'Wr 2 + 2rVz + V 2 ) d _ x , 
(z) = / r l — rVr a L dr 

Jo Vr 2 + 2rVz + V 2 

is negative and monotone non-decreasing; in particular, 

as z > 0, g'(z) > g'(-z). (A. 3) 

Using that 

T~ x u = T-^u = (cos^Wi + smcpu) d , w 2 ,---,Ud-i> — sin^Wi + cos ipuj d ), 
from l|A.l|) and (|A.2(1 one obtains 



p + (ta,mp) = / g(~ sin (puj 1 + cos (p to d ) dH d (A. 4) 

d— 1 



Now, substitute "— " for e and consider the orthogonal reflection U v with 
respect to the hyperplane {sin | wi = cos i| us d }; for any ^ G R d one has 

= (cosi/j^i + sint/j^, f 2 ,...,J/(i_i, sin ip vi - cos tp v d ) . 
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Changing the variable U^v = oj, one gets 

p_(tan^) = - / u 2 piU- 1 ^) dn d -\Lo). 

Using that U^ 1 = U v , one obtains 

p_(tany?) = — u)\ g(sintpuji — cosipujd) ctH (u>). (A. 5) 

The formulas l|A.4|) and (|A.5(1 can be written in the unified form 

p e (tamp) = e / lo 2 , g(e(— sin^^i + cos(pu>d)) dTi d ^ 1 (uj). (A. 6) 

Jst- 1 

Substituting ip = tt/2 in l|A.6|) . one obtains 

lim p e (u) = Pc(+oo) = e I uj 2 d g{-eu>x) dTi 4 " 1 ^). 

u — >+oo J S d ~~ 1 

Using that is invariant with respect to reflection lo\ h- > —u>i, one concludes 
that p + (+oo) = — p_(+oo), so (a) is proved. 

Further, using l|A.6|l . one concludes that the function p e is continuously 
diffcrentiable, and 

f d 

Pg(tani^) = ecos 2 ip ■ / uj d — — g{e(— sintpuji + cospujd)) dTi ~ (w) = 



dp 

= — cos 2 <p ■ / u j (cos tpuji + sm ipt^d) g\s(— sm ip uji + cos ipuid)) dH^ 1 ^). 
Jst- 1 

(A.7) 

Substituting ip = in HA.7|I . one obtains 

p'M = - [ J(ew d ) dH d -\uo). 

Js d s x 

Using that S^ -1 is invariant and the integrand is antisymmetric with respect to 
reflection u>i i— > one concludes that p' £ (0) = 0. Next, substituting ip = n/2 
in (|A.7(1 . one obtains 



lim p' £ (u) = p' £ (+oo) = 0- / i l4q'{-eu}{) dH d_1 (w) = 0. 



SI 

Thus, (b) is proved. 
Further, one has 

p' + (ta,np) -p'_(taxnp) = cos 2 p ■ / u; 2 d — $(tp, w d ) dH d ' x (u), (A.. 

Js d ~ 1 
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f d 

p' + (ta,n tp) = cos 2 tp- ujj —^> + (tp,uj 1 ,uj d ) dn^icj), (A.9) 

J S d-i dip 

where 

&(tp, u>i,C0d) = q{— sin tp u>\ + cos tp ujd) + g(sm tp u>i — cos tp u>d)- 
<&+(tp, wi,Wd) = g(- sin tp uji + cos tp uo d ). 

Designate 

f 9 d 

Jr c dtp 

!+ (c, tp) = UJ d — $+ (tp, uji , uj d ) dH 1 (uji , uj d ) , 

where r c = {(c^i, u>d) '■ uj2 + UJ d = c2 j < 0}- Let us prove that 

for any c e (0, 1) and tp £ (0, tt/2), Z(c, <p) < and I + (c, tp) < 0; (A.10) 

then, integrating 1(1 — |cl>| 2 , tp) and — p| 2 , y>) over a) = (u>2, . ■ . , Wd-i) and 
multiplying by cos 2 tp, one will conclude that the right hand sides of l|A.8|) and 
of l|A.9|) are negative, and so, (c) and the first inequality in (d) are true. 
Parametrize the curve T c by u>i = ccos9, Lj d = — csin0, 9 £ [0, 7r], then 



I(c, tp) 
where 



rir q 

/ c 2 sin 2 6» — \g(csin(tp + 9)) + g(-csm(tp + 6))] cd9 = <?T X +c 3 J 2 , 
Jo dtp 

Ji = / sin 2 <9 — L..1 d0, (A.ll) 

Jo dtp 1 ■> 

J 2 = / sin 2 ^ [...] d0, (A.12) 



and 
where 



/ c 2 sin 2 9 — g(-csin(tp + 9))cd9 = c 3 l+ + c 3 ^ 
Jo dtp 





p7T — 2ip 








Jo 






n = 


f ■ 




J 7T — 2ip 



(A.13) 
(A.14) 



Changing the variable ip = 9 + tp — 7r/2in IjA.f ljl . one obtains 

1\= \ cos 2 (tp — ip) — \g(c cosip) + g(—c cos ip)] dtp, 

J-tt/2+ v dy 



36 



and using the fact that the function ■ ■} under the sign of integral is odd, 
one gets 

I\ = / (cos 2 (ip — tp) — cos 2 (<p + tp)) -77 [g (ccostp) + g(—ccosip)] dtp. 
Jo dtp 

One has cos 2 (y — tp) — cos 2 (if + tp) = sin2<ys sin 2-0 > 0. Taking into account 
HA.3|) . one also has that i [•••] = — csinV> (g'(ccostp) — g'(— ccostp)) < 0. 
Hence, I\ < 0. 

Making the same change of variable in l|A.13|) , one gets 

It = / cos 2 (ip — tp) — g(— ccostp) dtp = 

J-ir/2+ V dtp 
ir/2-cp j 

(cos 2 (f — tp) — cos 2 (ip + tp)) — g(— ccostp) dtp. 

dtp 

One has cos 2 (tp—tp)— cos 2 (ip+tp) > 0, and g(— ccostp) = c sin tp g' (—ccostp)) < 

0, thus 1+ < 0. 

On the other hand, changing the variable x = 9 + tp — 7r in (|A.12fl . one 
obtains 

I2 = J sin 2 (> - x) ^ [£»(csinx) + £>(-csinx)] dx- 
The function 3^ [• • •] is °dd, therefore 

(sin 2 (ip - x) - sin 2 (tp + x)) ^- [p(csinx) + g(-csinx)] c?X- 

One has sin 2 (</? — x) ~ sin 2 (tp + x) — —sin 2ip sin2x < 0, and -j- [...] = 
ccosx (g'(csinx) — q'(— csinx)) > 0. Hence, I2 < 0. 

Further, as ir — 2ip < 9 < it, one has g(— csin(<y9 + 9)) — —ccos(ip + 

9) g' (—csin(ip + 9)) < 0, and using (|A.14Jl . one concludes that 2"^ < 0. 

Thus, the inequalities in IjA.lOfl are proved, and so, (c) and the first inequality 
in (d) are true. 

Passing to the limit ip — > tt/2 in l|A.5|l . one gets 

d-l 



lim p_(tan<p) = p_(+oo) = - / uj 2 d q(ui) dH (uj). 

If—- >ir/2 



Thus, to prove the second inequality in (d), one needs to verify that 

ui\ g(u>i) dH^ 1 ^) > I uj 2 d ^(sin^cji - cosLpuj d ) dH d ^ 1 (uj). 

d-l J gd-l 

Denote 

J(c,ip) = / i^ 2 g(sin ip L1J1— cos (p ujd) dTi 1 (uJi,ujd) — c 3 / sin 2 9 g(csin(ip+9)) d9. 
Jr c Jo 
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It suffices to prove that 

for any c G (0, I) and <p G (0, tt/2), J(c, tt/2) > J(c, <p); (A.15) 

then by integrating J(l — |o)| 2 ,7r/2) and J(l — |w| 2 ,(^) over £> = (u 2 , ■ ■ ■ ,^d-i), 
the inequality (|A. 15|> will be established. 

One has J(c, (p) = J\ + J 2 , J(c, n/2) = J-j" + J|, where 

/•TT/2 — if />7T 

Ji = / sin 2 6»£>(c sin(</3 + 0)) d6>, J 2 = / sin 2 6» g(csin(( / 3 + 6>)) d9, 

/•ir/2+ip p-it 

J{= I sin 2 9 gic cos 9) d9, J 2 = I sin 2 9 g(c cos 9) d9. 

JO Jir/2+ip 

As < 9 < 7r/2 — if, one has — cos# < < sin(p + 9), hence g(—ccos9) > 
g(csin(ip + 9), thus 

p-n/2 — ip pir/2 — 

J 2 = sin 2 6> g(-ccos6l) d# > / sin 2 6> p(csin((p + 9)) d9 = J x . 

Jo Jo 

(A.16) 

Further, one has 

i-Tr/2+ip i-Tr/2+ip 

2J{ = sin 2 9 g(c cos 9) d9+ sin 2 (n/2+tp-9) g(ccos{n/2+(p-9) d9, 

Jo Jo 

t-ir/2+ip rir/2+ip 

2J 2 = sin 2 9 g{csin{9-tp)) d9+ sin 2 {ir /2+>p-9) g{csin{ii / 2-9)) d9, 

Jo Jo 

hence 

r-K/2+ip 

2J 1 *-2J 2 = / [sin 2 6- cos 2 (9 - <p)][g(c cos 9) - g(csm(6 - <p)]d6. (A.17) 

Taking into account that the function g is monotone decreasing and that sin 2 9 — 
cos 2 (9 — (p) — (sin(6> — ip) — cos6 l )(sin(0 — <p) + cos 9), one concludes that the 
integrand in l|A.17|) is positive, hence J I > J 2 . From here and from l|A.16|) it 
follows that l|A.15|) is true. Lemma 1 is completely proved. □ 

Proof of lemma [5] 

(a) Parametrize the set S_ according to v\ = — sin6>, v 2 = —cos 9, 9 G 
[-tt/2, tt/2], then HA.6|) takes the form 

p £ (tantp) = e cos 2 9 g(-ecos(ip + 9)) d9, e G {-, +}. (A. 18) 

J -tt/2 

Twice differentiating both parts of this equation with respect to <p>, one obtains 
^'(tan^) _ _ r' 2 d 



. , cos 9 2 sin ip — g(—ecos(ip + 9)) d9+ 
cos 6 ip J_ 7! /2 op 
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-£ 



rn/2 q2 

I cos 2 9 cosp — — ; - g{— ecos(p + 9)) dd. 

J-n/2 Oif 2 



t/2 

Integrating the second integral by parts and taking into account that -J=^ g(~s cos(tp+ 
0)) = Wk g(-ecos(tp + 6)), k = l, 2, one gets 

PeK - V) =e 2 cos 9 sin(0 - ip) — g(-e cos(p + 9)) d9. 
cos 3 ip J_ n/2 99 

Integrating by parts once more and denoting g E {jf) '■— p"(tan<p)/(2cos 3 ip"), one 
gets 

/■t/2 

g e ((p) = -s g(-s cos(cp + 6)) cos(26-ip)d6 . (A.19) 

Fix the sign "+" and prove that 

(I) for < ip < tt/6, g+{p) < 0; 

(II) for p > 0.3tt, g + {ip) > 0; 

(III) for n/6<p< 0.3tt, c^) > 0. 

The relations (I), (II), and (III) imply that there exists u+ G (l/\/3, tan(0.37r)) 
such that p'+{u) < as u € (0, m+), and p+(w) > as u £ (u + , +oo). 
(I) Changing the variable ip = — ip/2 + tt/4, one gets 

/ 371-/4-^/2 
2(-cos(3<p/2-7r/4 + V>)) sin 20 (A.20) 
-Tr/i-ip/2 



L\ + £2, where £1 



■n/A+ip/2 p3ir/A-ip/2 

■ ■ • , £ 2 = - / 

-ir/A-ip/2 Jn/A+ip/2 



One has 



/.tt/4+^/2 

£1 = / [g (- cos (3^/2 - tt/4 - ^)) - £ (- cos (3^/2 - tt/4 + 0))] sin 2^ #. 

Jo 

(A.21) 

One has < 2ip < ir/2 + ip < tt, hence sin 2^ > 0. Using that < ip < tt/6, 
one obtains that — tt/2 < 3p/2 — tt/4 — ip < — \3<p/2 — tt/4 + ip\, therefore 
— cos(3<ys/2 — 7r/4 — ip) > — cos(3ip/2 — tt/4 + ip). The function g monotone 
decreases, therefore g (— cos(3ip/2 — 7r/4 — -0)) < g (— cos(3</?/2 — 7r/4 + ^)). 
Thus, the integrand in HA.21|) is negative, and so, C\ < 0. 

Change the variable x = ip — 7r /2 in the integral £2- One obtains 

r-K/A-tp/2 

C 2 = g(- cos (3p/2 + tt/4 + X )) sin 2 X d X = 

J -Tt/A+ip/2 

= / [g (- cos {ip/2 + tt/4 + x)) - g(- cos (3p/2 + tt/4 - x))] sin 2 X d\- 

Jo 

(A.22) 
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One has < 3^/2+tt/4-x < 3ip/2+7r/4+x < n, hence - cos(3^/2+7r/4-x) < 
-cos(3v3/2 + 7r/4 + x), and so, q{- cos(3^/2 + tt/4-x)) > q{~ cos(3ip/2 + 7r/4 + 
x))- Therefore the integrand in l|A.22|) is negative, and £2 < 0. Thus, (I) is 
proved. 

(II) By lpO0)l . one has 

g+(<p) = Ii + 1 2 + + 1a, 

where 

/<p r—<p rTr/4+tp/2 />37r/4— ipj 

(...), l2 = ~ (...), 1 3 = - (...), lA = - 

-ip J -ir/i-ip/2 J tp Jtt/4+ip/2 

and (...) = q{— cos(3iy9/2 — 7r/4 + ?/;)) sin 2-0 dip- 
One has 

Xi = / [e(-cos(3^/2— 7r/4--0))-e(-cos(3(y9/2-7r/4+0))] sin2^d<0; (A.23) 
^0 

using that 93 > 0.37T, one easily verifies that as < ip < ip, \3ip/2 — 7r/4 — ^| < 
3^/2 - 7r/4 + t/» < 7r, hence - cos(3</?/2 - 7r/4 - ip) < - cos(3^/2 - 7r/4 + ^). 
Using that g monotone decreases, one concludes that the integrand in (|A.23ll is 
positive, thus X\ > 0. 

Next, as — 7r/4 — ip/2 < ip < —ip, one has sin 2ip < and |3<^/2 — 7r/4 + ip\ < 
iip/2 - 7r/4 + ip + 2(p < 7r, hence - cos(3<p/2 - 7r/4 + 0) < - cos(3<p/2 - 7r/4 + 
ip + 2(f), and thus, g(- cos(3(^/2 - vr/4 + ip)) < g(- cos(3^/2 - tt/4 + + 2ip)). 
Therefore, 

p(- cos(3<^/2 - tt/4 + V + 2ip)) sin 20 d0 = 

-7T/4-V/2 



x/4+<p/2 

g(-cos(3^/2-7r/4-x + 2^)) sin2xdx, 



and 



p?r/4+y/2 

l2+?3> I [g(-cos(3<p/2-ir/4-ip + 2tp))- 



-g{~cos(3ip/2- n/4 + ip))} sin 2-0 #. (A.24) 
On the other hand, one has 

<.3Tr/4-^/2 

X 4 = / [e(-cos(3^/2+37r/4-V'))-e(-cos(3^/2-7r/4+0))] sin2ipdip. 

Jtt/2 

(A.25) 

Changing the variable 9 = ip — ip in l|A.24|) and 9 = ip — tt/2 in (|A.25|) and 
summing both parts of these relations, one obtains 

7r/4-y/2 



1 2 +1 3 +1 4 > / ^{9)d9 
Jo 
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where 

= [Q (- cos(5<^/2 - tt/4 - 0)) - g (- cos(5^/2 - tt/4 + 0))] sin(20 + 2tp)- 

-[£»(- cos(3(^/2 + 7r/4- 0)) - g (- cos(3^/2 + tt/4 + 0))] sin 26*. 

Let us show that *(0) > 0; it will follow that X 2 +^3 +Z* > 0, and thus, (II) 
will be proved. 

One has < 26* < tt/2 - tt/2 - ip < 2ip < 29 + 2<p < tt/2 + ip, hence 

< sin 20 < sin(20 + 2(f), (A.26) 

Denote J x (0) = g (- cos(5<^/2 - tt/4 - 0)) - g (- cos(5<^/2 - tt/4 + 0)), J 2 (0) = 
£> (- cos(3<^/2 + tt/4 - 0)) - £ (- cos(3^/2 + tt/4 + 0)). One has 

Ji(0) = - [ Q'{- cos(5^/2 - tt/4 + X )) sin(5^/2 - tt/4 + X ) d X , 
^(0) = - / g'(- cos(3<^/2 + tt/4 + x)) sin(3^/2 + tt/4 + X ) d X ■ 



Using that ip > 0.3tt, one gets ir — 2ip < 3ip — tt/2 < 5ip/2 — tt/4 + X < 2<p and 
2ip < 3(f/2 + tt/4 + X <p + tt/2 < tt, hence 

sin(5<^/2 - tt/4 + X ) > sin(3v3/2 + vr/4 + X ) > (A.27) 

and — cos(5iy9/2 — 7r/4 + x) < — cos(3<^/2 + tt/4 + X ). Taking into account that 
g' is negative and monotone increasing, one gets 

g'(- cos(5(y5/2 - tt/4 + X j) < g'(- cos(3<p/2 + tt/4 + X )) < 0. (A.28) 

From HA.27|) and (|A.28|I it follows that J x {9) > J 2 (9) > 0, and taking into 
account (|A.26|) . one concludes that , 3>(0) > 0. 
(Ill) One has 

g' + {tp) = --T- / £>(-cos(3</2/2 - 71-/4 + ip)) sin2ipdip = 

UP J-- K /i-ip/2 

3tt/4- V /2 q 

— g(-cos(3ip/2-Tr/4 + ip)) sm2ipdip+ 

-Tr/4-^/2 Oip 

+ — cos ip [q{— cos(ip — tt/2)) — g(— cos((p + tt/2))]. 

One has cost/? > 0, — cos(p — tt/2) < < — cos(p + tt/2), hence g(— cos(<p — 
tt/2)) - g{— cos(p + tt/2)) > 0, and thus, 

/3tt/4-¥>/2 Q 
— e(-cos(3<p/2-7r/4+?A)) sin2^# = -(£i+£ 2 +£ 3 +AC 4 ), 
-7r/4-<f./2 dL P 
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where 

/n/2-3ip M ,-ir/A+ip/2 ,-3v/i-ip/2 

K.2 = / /C3 = / IC4 = 
-n/A-tp/2 Jir/2-3ip JO Jn/4+ip/2 

In order to prove (III), it suffices to show that JC\ < 0, JC2 < 0, /C3 < 0, and 
JC 4 < 0. 
One has 

3 r-Tv/2 — 3tp 

JCx = - / cos(3<^/2 - tt/4 + V)) sin(3^/2 - tt/4 + ijj) sin 20> #. 

2 J-Tr/4-<p/2 

Using that (p < 0.37T, one easily verifies that — 7r/4 — ip/2 < n/2 — 3ip. For 
-tt/4-^/2 < -0 < 7r/2-3i^, one has -tt/2 < 3<p/2- tt/4 + < 0, -tt < 20 < 0, 
hence sin(3i/?/2 — 7r/4 + -0) < 0, sin 20 < 0, and taking into account that g' < 0, 
one concludes that JCi < 0. 

Changing the variable = 3tp/2 — 7r/4 + 0, one has 

3 r3 V /2-w/4 d 

JC 2 = - / — p(-cosx) cos(2 X - 3tp) dx = 

Z Jn/4-3ip/2 "X 




— g{- cos x) [cos(2x - 'iip) - cos(2x + 3<p)} d\- 



One has ^ g(— cos x) < 0, cos(2x — 3<p) — cos(2x + 3<p) = 2 sin 2\ sin 3y > 0, 
hence Ki < 0. 

Further, one has 

JC 3 = - / e'(-cos(3(^/2 - tt/4 + 0)) sin(3^/2- tt/4 + 0) sin2V>#. 

2 Jo 

It is easy to verify that the integrand is negative, hence K.3 < 0. 

Changing the variable = — 7r/2in the integral JC4, one obtains 

[Tt/i-ip/2 g 

JC 4 = - — g(— cos(3(^/2 + tt/4 + (9)) sin 20 d0 = 

J-w/4+<p/2 Of 




■K/i-tp/2 

[g'(- cos(3<^/2 + tt/4 - 0)) sin(3^/2 + tt/4 - 6)- 



-g'(- cos(3</?/2 + tt/4 + 0)) sin(3<p/2 + tt/4 + 0)] sin 29 d9. 

Using that ip > 7r/6, one easily verifies that {if/2 + n/4 + 9) — tt/2 > \(3(p/2 + 
7r/4-6»)-7r/2|, hence sin(3< y 9/2+7r/4+6») < sin (3^/2+ tt/4 -0) and - cos(3^/2 + 
tt/4 + 6) > - cos(3(^/2 + tt/4 - 9), therefore < -g'{- cos(3<^/2 + tt/4 + 9)) < 
-g'{- cos(3^/2 + tt/4 - 9)). This implies that /C 4 < 0. 

(b) Let us slightly change the notation just introduced: we shall write down 

/>oo 

g(z, V) = r 3 o{^Jr 2 + 2rVz + V 2 ) dr, 
Jo 
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Figure 11: 



/■tt/2 

p_(tan<y9, V) = - cos 2 9 g(cos(tp + 9), V) d9, 

J-tt/2 

g-{(f,V) = - rT J„ = tan V P-(u,^)/(2cos 3 (^) = / g(cos(<p+0),V) cos(20-(p)d9, 



du 



tt/2 



thus explicitly indicating dependence of these functions on V. Using these 
formulas, we see that the function a a ^{r) = a{r) + aa((3r) generates the 
function g a '^(z) = g(z,V) + -^g(z,/3V) and the pressure function p_(u) ~ 
p_(u,V) + $p-(u,f3V). 

Suppose that the set O := {u : p-(u, V) > p-(u, V)} coincides with an 
interval (0, vP_ ); otherwise, the hypothesis of lemma|21(b) is valid for a = and 
arbitrary (3 > 0. Consider two arbitrary values u\ S (0, u_) and it 2 > vP_ , and 
designate 

p_( Mi ,vO-p_(o,vq p_( M2 ,y)-p_(o,y) 

Ai(p) = , A 2 [p ) = , 

Ul u 2 

Ai 2 (p) = ; 

u 2 - Ul 

one has A x (p) > A 2 (p) > A i2 (j?) and p'_(tt 2 , V) > A 2 (p) (see Fig.H]}- 

Note that the function g~(<p, V) is continuous with respect to <p on [0, tt/2), 
and has a limit as <p — » 7r/2 — 0, hence the value g := sup ¥3g [ ,7r/2) \9-(<Pi ^01 i s 
finite. Denote 



u = mm 



Aip-A i2 p p'_(u 2 ,V) - A 2 p 



2 ' 2 

It is easy to see that if the functions p(u) :— -^p-(u, (3V) and g((p) := 
ft -£fi\u=t a n V P-{u,PV) satisfy the inequalities 

\p'(u)\ < lu for < u < u 2 , p'(u) > -CJ for u > u 2 , (A.29) 
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and 

g(ip) < — g for some <p > arctanw 2 , (A. 30) 

then the function ptl^ \u) — p~(u) + p(u) has the following property: the set 
Qa,p _ | M . <p°^P(y^ > p" ,,3 (u)} has at least two connected components; one of 
them is contained in (0, u 2 ), and the second one contains tantp and is contained 
in (u 2 , +00). 
Denote 

P fi (ip) =mscx{ sup \p'_(u,/3V)\, sup (-p'_(u, (3V))}. (A.31) 

0<u<tanip u>tanip 

The rest of the text (till the end of Appendix A) is devoted to the proof of the 
fact that for any s > there exist <fr, ip, and /3 such that 

vr/2 - e < < vr/2 and < ——7—— < e; (A.32) 

then, letting <j> = tan«2, £ = u/g, one can find a such that l|A.29() and (|A.30|I 
are valid, so the statement of lemma 121(b) for given a and (3 is fulfilled. 

Let us carry out some auxiliary calculation. For <p <G (7r/10, 7t/2), (3 > 
define R(<p, (3) by 

Then pick out constants a > 1 , c > in such a way that there exists a circular 
sector £ of angle c, outer radius R, and inner radius i?/a, with the center 
(0,—/3V), that belongs to the angle U :— {v = (di,^) : «i > 0, — tan(^ — 
-£) < f 1 < — tan(| — <£>)} (see Fig.^J. One can choose, in particular, arbitrary 

ce (0, tt/5) and a = c c ™ff ■ 

Next, using the definition of g and passing to the variables i>i = rsinx, 
v 2 = rcosx, one obtains 

<p+i r r 

g(cosx,l3V)d X = (vl+vl)a{\v + (3Ve 2 \)dv 1 dv 2 . (A.34) 
^+2= JJu 

Taking into account that for v € £ one has u 2 + u| > — ^) 2 > — 
Rj sin(-^ L — and S CM, and passing to polar coordinates, one gets 

(v\ + v%) a(\v + pVezl) dvidvi > 



>R 2 . ^ -1 <y(\v + PVe 2 \)d Vl dv 2 



iS in 2 (f -f) 



1 




Lin 2 (f-f) J V 
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Figure 12: 
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> ( ■ 2(5 l 5n - 1 J ' c I* r 3 a(r) dr. (A.35) 

From <|A.34|) and l|A.35|) one concludes that 

rv+f f 1 \ rR 

g(cos X , PV) d X > g— - 1 • c / r 3 a(r) dr. (A.36) 

\ Sm ll 4 ) / 

Besides, one has 

^ 4 / 4 ^ 8 

^(cos x, /3V) dx<2 g(cos x, 0V) dx < 

< fj {vl +v%)a{\v+i3Ve 2 \)dv 1 dv2, (A.37) 

where V = = (^1,^2) : \v + f3Ve 2 \ > R} is the complement of the circle 
of radius R with the center (0, — j3V). Taking into account that for v G V, 
0V < \v + /3Ve 2 \/sm(^-^), one has |«| < pV + \v + 0Ve 2 \ < \v + f3Ve 2 \ (1 + 
l/sin(f-^)), hence 



(vj + vj) a(\v + f3Ve 2 \) d Vl dv 2 < 



2 

V 

\ 2 



- ( sin(^-5se) +1 ) jJ^ + PV^l °(\ v + PVe 2 \)d Vl dv 2 



Vsin(f - 
From iA~37l) and llXT38l) one obtains 



\ 2 00 

t — + l) -2% r 2 a(r) rdr. (A.38) 

T - 4 ) / •>« 



2^4 / I 



p(cos x, /3F) d X < I — ^ g— + 1 ■ 2tt / r J a(r) dr. (A.39) 

k sin(-g f) J Jr 

Further, by hypothesis (b) of lemma [21 f° r n > and for r sufficiently large 
the function 7(7") := r n+3 a(r) monotonically decreases, hence for R sufficiently 
large one has 

1 n—n+l rR R~ n+1 

r~ n 1 (r)dr < 1 {R) -. / r~ n 7 (r) dr > 7 (i?)(a"~ 1 - 1) , 

r n - 1 J R/a n - 1 

thus 

/•oo 2 r^- 

r 3 tr(r) dr < : / r 3 cr(r) dr. 

r a™ 1 - 1 J R/a 
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By virtue of arbitrariness of n, one concludes that 



r 3 a(r) dr 

lim — = 0. A.40 



R^+oo J" r 3 cr(r) dr 



From (fP0)l . (fQ6)l . (fQ9l) . and l|A~33)> it follows that for any tp e (tt/10, tt/2) 



Lli 4 e(cos X ,/3y)d x 
lim ^+1 — ; = o (A.41) 



and 



g(cosx,0V)dx>Ci(if) / r 3 cr(r)dr, (A.42) 

^ + ^f J R/" 

where ci(v?) = c(l/sin 2 (^-^)-l). 

Substituting "— " for e and changing the variable x — f + $ m (IA. 191) . one 

gets 

5 _ fo>, /3V) = / g(cos X , /3V) cos(2 x - Zip) d X - (A A3) 

Jip—ir/2 

Denote U\ := -j> = (vi,v 2 ) : jf^- > -tan<^}, Vi := -j> = (vi,v 2 ) : |v + /3Ve 2 | > 
/3 V cos yj}; Vi is the complement of the circle of radius (3V cos </? with the center 
(0,-PV). One has Ui C Vi, and for v G Vi |u| < |w + /3Ve 2 | + /3V" < |<u + 
/3 V e 2 1(1 + 1/ cos y>) , hence 

|5-(VJ,j9V)| < r +V %(cos X ,/3y)d x - // (« 2 +« 2 V(k + /3Ve 2 |)ch; lC fc 2 < 

J-ip-Tr/2 J JUi 

1 ^ 



<(1 + ^^J yy |« + /3Ve a r tr(|t; + j9Ve 2 |) dvidv2 



1 ^ 2 



oo 

2 



1 H • 2tt / r 2 cr(r) rdr. (A.44) 

COSy)/ Jpvcosip 

Taking into account that cos(2x - 3<ys) < as x e [if + f , x + cos(2x - 
3y) < cos -|) = - cos (| + f) < as x e [^ + ^-,ip+ §], and using 
<|A.43(1 . one obtains 



2^4 



{<p,0V)< ^ e(coBX,/9V r )dx-cos(J + |). / R Q(axx,0V)d X . 



4 ~ 8 



(A.45) 

From lf06| . (S33J, l|A~42)) - and <fA~45|) it follows that for any y and for /3 

sufficiently large 



R 

3 



9-{<p,PV)<-o»{<p) / r A a{r)dr, (A.46) 
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where R = R(<p,(3) and C2{<f) = \ cos (f + ' c i(v)' 

Next, con sider = 4>((p) := - §; one has ^Fcos^ = /3Fsin + ^ ) 
i?; hence, by (|A.44jl . for any ip e [0, 4>] 



1 ^ 2 



cos ip 



\g-(4>,0V)\< 1+ 2tt / r^(r) dr < 











\ COS0 





/3V^ cos tp 

7{r) dr. (A.47) 



Further, one has 

/■arctan u 



p'_ (u, 13V) = / p" (v, (3V) dv = 2 / 3 _ (V», /?V) cos ^ 

Jo Jo 

hence, according to (|A.47|1 . 

sup \p'_{u,f3V)\ <2 / \g_y,,0V)\dil>< 

O<«<tan0 JO 

< 2(f) ( 1 + —^—) 2tt f r 3 a(r) dr (A.48) 

V COS0/ J R 

and 

sup (-p'_(u, f3V)) = - / j£(i/,/3V)dH- sup I- / p'^(i/,/3F)c2z/) < 

u>tan0 JO -it > tan <f> \ Jtan^ / 





1 












COS 





r 3 <r(r) dr + 2(tt/2 - 0) ■ sup 

V>6[0,tt/2] 

(A.49) 

Using (|A.31fl . (|A.48|) . HA.49|) . and (IA.46f) . one gets that for [3 sufficiently large 

0< ^ <C3M |^ + »/2-ft (A.50) 

sup^ e[0i7r/2 ] (-.9- W>, pV)) J« r 3 CT ( r ) dr 

where i? = #(<£,/3) and c 3 ((p) := </> (l + ^) • For arbitrary e > 0, 

choose y such that < 2(ir/2 — <fr(<p)) < e/2, and then, taking account of 
<|A.40(1 . choose (3 such that the first term in the right hand side of (|A.50|I is less 
that e/2. The relation (|A.32|) is established. □ 

Appendix B 

d = 2 Here the formula H2.10[l takes the form 



(viu + ev2)3 
1 + u 2 



p e (u, V) = e I I TT~2 P( v ) dvidv 2 , e e {-, +}. 
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Passing to the polar coordinates v = (— r sin ip, —er cos <p) and using that p(v) 
a(\v + Ve2\) — <j(r) — eV cos <p a'{r) + o(V), V — > + , one obtains 



o y i + u 2 



p £ (w, V) = e J J — ; — 2 ((r(r)-eV cos <p a {r)+o(V)) r dr dip 



a(r) r 3 dr ■ 1^ -V a'(r) r 3 dr ■ + o(V), (B.l) 



where 

1^ = / (cos(p — p>o))+ dip, = / (cos(ip — <po))+ 2 cos ip dip, 

Jo Jo 

x + := max{x, 0}, ipo '■= arccos -/jjcp ■ Changing the variable j/j = 93 — </Jo ii 
these integrals, one obtains 

rir/2 A 

J (2) = / cos 2 ip dip = 7r/2, J (2) = cos </? / cos 3 i/"# = 



-7T/2 J-tt/2 3V1 + M 2 

Substituting the obtained values in l|B.l|) and using that — J Q °° cr'(r) r 3 d?' = 
3 J °° cr(r) r 2 dr, one comes to the formula (|4.3.2|) with coefficients Q4.8[l . 

d = 3 Formula l|2.1U[) takes the form 

p e (u, V) = e J J J ^-j-q^^— p{v)dv!dv2dv 3 , e G {-, +}. 

Passing to the spherical coordinates w = (— r sin</?cos#, — r sin (/j sin 9, —er cos (/?), 
one obtains 

/ t r\ [ 27T f 00 r 2 (shop cos9u + cos p) + \ , ^ 
p e {u,V)=e / / — — ^- (a(r) - eV cosipa'(r)+ 



a Jo Jo l + u 



+o(V))r 2 drd9 sin ipdtp = e \ a(r)r 4 dr-I® -V \ a 1 (r) r 4 dr ■ J (3) + o(V) , 

(B.2) 



where 



r(31 /-^r 2 1 + U 2 COS 2 . Ja 

= / / (cos(^-0o(^)))+ — — -5 — simpdipdd, 
Jo Jo l + u 

r(3) / 2l r, , - 2 1 + U 2 COS 2 . , Ja 

/ / (cosC^ — ( Po(c')))+ 2 cos ip simpdipdd, 

Jo Jo l + u 

<po(0) = arccos V / 1+M 2 cos 2 g ■ Changing the variable -0 = 93 — <Po(@)> °ne obtains 

/(3) = / 1 + " C ° 9 b sin O (g) / cos 3 ^ = 

Jo 1 + u J -Me) 
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27r 1 + u 2 cos 2 6> (1 + sin0 o (0)) 2 



1 + u 2 3 
271 (y/1 + u 2 cos 2 9 + u cos e f 



dO = 



o 



3(1 + u 2 ) 



j( 3 ) = / de 



2tt 



1+u COS 
1 + u 2 



2 nr> c 2 a rir/2 



,„ 1 + u 2 cos 2 f 
o 2(1 + u 2 ) 

tt/2 



2 sin(2V> + 20 o (fl)) 
cos V ^ «V J 

cos20o(#) / cos 2 -0 sin 2ip dip+ 
J -Me) 



■sin 200(f) 
1 



cos tp cos 2-0 c£0 



2k 



2tt 



1 + u0 o (6>) cos 9 



i(l + u 2 )J 



o 4(1 + u 2 ) 
[1 + (itcos#) arctan(ucos#)] dO 



d9 



2 % /TT^ 2 



Substituting these values in (|B.2|I and using that — J °° cr'(r) r 4 dr — A J °° <r(r) r 3 dr, 
one gets the formula (|4.3.2|l with coefficients (|4.9() . 
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